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A STUDY ON CHARACTERIZATION OF STRATOSPHERIC 
AEROSOL AND GAS PARAMETERS WITH THE SPACECRAFT 
SOLAR OCCULTATION EXPERIMENT 

By 

William P. Chu 1 
SUMMARY 

This report contains results of research investigations performed during 
the period from October 1975 to December 1976 under NASA Grant N3G 1259, 
Spacecraft remote sensing of stratospheric aerosol and ozone vertical profiles 
using the solar occultation experiment has been analyzed. A computer 
algorithm has been developed in which a two-step inversion of the simulated 
data can be performed. The radiometric data are first inverted into a vertical 
extinction profile using a linear inversion algorithm. Then the multiwav»'e- 
length extinction profiles are solved with a nonlinear least-square algorithm 
to produce aerosol and ozone vertical profiles. Examples of inversion results 
are shown illustrating the resolution and noise sensitivity of the inversion 
algorithms . 


INTRODUCTION 


The minor constituents in the stratosphere have received great attention 
in recent years because of concerns over their roles in affecting the global 
environment. With mankind's increasing activities in the stratosphere such as 
SST and Shuttle missions, and ground-based activities such as those involved 
with nitrogen fertilizers, an increase in stratospheric aerosol loading 
together with the removal of stratospheric ozone are possible through various 
chemical reactions. A realistic assessment of this problem is still not 
possible due to the lack of information on. the distributions of most of the 
constituents in our atmosphere . 
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This report analyzes one of the atmosphen ' remote sensing methods using 
spacecraft solar extinction technique to retrieve stratospheric aerosol and 
ozone vertical profiles. In this remote sensing technique, the spacecraft 
instrument will direct towards the horizon and measure transmitted solar 
radiant intensity in a number of spectral intervals from 0.38 to 1.00 micron 
wavelength. The locations of the spectral channels are selected based on 
considerations of maximum extinction contributed from any individual consti- 
tuent, interference from other constituents, and available energy for 
measurements. The radiant intensity data are then inverted to produce vertical 
profiles for each of the constituents. The inversion procedure can be 
separated into two steps. The first step involved retrieval of the total 
vertical extinction profiles at each spectral interval. The second step then 
separated the contributions from each constituent using the multispectral 
extinction data at each altitude level. Examples of inverted results from 
simulated data including experimental errors are presented and discussed in 
this paper. 


EXPERIMENTAL CONCEPT 

In the solar extinction experiment, multi wave length measurements of the 
transmitted solar radiant intensity along the line of sight through the 
atmosphere are made on board the spacecraft. The basic geometry is illustrated 
in figure 1. As the spacecraft emerges from Earth's shadow, the radiometer 
aboard the spacecraft will point to the sun, measuring solar radiant intensity 
at several preselected wavelength regions. In this paper we will consider a 
radiometer with a small field of view in order to achieve higher vertical 
resolution. The radiometer is assumed to either lock on a fixed position on 
the solar disk or scan across the solar disk for data acquisition. 

The monochromatic radiant intensity reaching the radiometer is given by 
the fundamental equation of radiative transfer (Chandrasekhar, 1950). 

W - W WV + / „ „ , •V h > d W h > CD 


L 


RADIOMETER 





Figure 1. Geometry of tne solar extinction experiment. 


The above equation relates at wavelength \ , the radiant intensity 
I^Ch^) measured at distance h 2 , to the radiant intensity I^(h o ) at 
boundary h Q , the monochromatic transmittance T^(h 2> h Q ) and source 
function J^(h) at distance h along path from h z to h Q . 

For the solar extinction experiment with spectral channels closed to the 
visible region, the contribution by atmospheric emission to the radiant 
intensity is very small and can safely be neglected. The radiative transfer 
equation is thus greatly simplified to 

W ■ W WV m 

where the radiant intensity I^(h o ) is now the solar radiant intensity, 
and T^(h z ,h o ) is the transmittance of the atmosphere between the sun and 
the spacecraft. 

In the wavelength range from 0,3 to 1.0 micron, atmospheric attenuation 
are predominantly caused by aerosol, ozone, and air molecules [Rayleigh 
component) with minor contributions from N0 2 , water vapor and 0 2 molecules. 
The model of the stratospheric extinction versus wavelength profile can be 
constructed from different vertical distributions of aerosol and ozone. 

Figure 2 shows an extinction coefficient model as a function of wavelength 
at an altitude of 18 km. Notice that at most of the wavelength shown, the 
total extinction consists of contributions from aerosol, ozone, and Rayleigh 
component with different weights. The weightings will also be different for 
different altitude levels, depending on the distributions of the different 
constitutents. 

If we made the assumptions that the atmosphere is spherically stratified 
and spherically symmetric, the transmittance function t. (h ,hj will then 
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be given by the Lambert- Beer law as 
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Figure 2. Atmospheric extinction versus wavelongth model at an altitud.* of 18 km. 


where 0^(h) is the atmospheric extinction coefficient at wavelength X 
versus vertical height h , P(h,h T ) is the optical path at height h for 
line of sight with tangent height h T . Equivalently, we can define the 
optical depth g^Ch^) as 


g^(h T ) = in U/T v (h z ,h o )} 

co 

= 2 / 0 x (h) dP(h,h T ) 
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The determination of vertical extinction proille at each wavelength 
X required the inversion of the integral in equation (4) . 

INVERSION TECHNIQUES 

Inversion of Measured Transmittance to Total Extinction Profiles 

In equation (4), the measured parameter is equal to a kind of 

convolution of the weighting function 3P(h,h T )/9h with the unknown vertical 
extinction profile 3^00 • The integral can be approximated with a sum 
over n discrete atmospheric layers with equal thickness and assigning to 
each layer an averaged extinction coefficient 0^ . The integral equation 
can then be replaced by a system of linear equations. 
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where P^. is the optical path length in the jth layer with tangent height 
at the ith layer. Equation (5) can be abbreviated into matrix form 


g = P3 


(6) 


6 


where g is the column vector for the measured optical depth, 0 is the 
column vector for the unknown extinction profile, and P is the optical 
path length matrix. 

Equation (6) can be solved directly for 0 either exactly when n » m 
or in the least square sense when n < m 


0 « P’ 1 g 

n = m 

(7a) 

0 = (pVr 1 p t g 

n < m 

(7b) 


T 

where P" 1 is the inverse of P and P is the transpose of P . The 
solutions obtained from the direct inversion are generally unstable due to the 
presence of noise associated with the measured parameter g . In order to 
suppress the unphysical oscillations in the inverted solutions, some con- 
straints on the high frequency components of the solutions have to be 

* 

incorporated in the inversion schemes. We have analyzed two different 
inversion methods. The first method is the linear constrained inversion as 
developed by Twomey (1963) and Philips (1962) ; the second method is the 
iterative scheme recently developed by Twomey (1975). We will discuss the 
two methods separately below. 

Linear constraint inversion 

The solution to equation (6) is given by the following expression 

0 = (P T P + yH)" 1 P T g (8) 

where normally H is the constrained matrix, and y is a constant smoothing 
parameter. In our analysis, we used the constrained matrix H which minimizes 
the second difference of the solution. The smoothing parameter Y whose 
magnitude is proportional to the noise level of the measurements has been 
replaced by a diagonal matrix of element Y^ , The reason for using a variable 
smoothing parameter Y arises from the nonlinear expression for the noise 
term in equation (4) . Assuming the radiometer has a noise level e^ 
associated with each radiant measurement 1^ . The deduced transmittance 
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will have approximately the same noise level e^ , By computing the optical 
depth , we have 


g i = In (1 / (t^ + e i )} 


C9) 


For e. << t. , we can expand (9) and obtain 
1 i 


“ In Cl/T i } - e^/r^ 
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The second term on the right side of the equivalent noise level for the 
parameter g. . Thus we choese the matrix element y. . as 
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( 11 ) 


where y Q is a constant parameter to be adjusted for optimum inverted 
solutions. 


Iterative inversion 

The iterative method starts by initially assumed a guessed solution, 
then the solution is continuously updated through the following iteration, 


pCk+i) = gCk) + ,^00 _ , gCh) 

P 3 3 i J 13 3 


( 12 ) 


where gP*) is the solution 0^ at k-iteration, yP^ is the ratio of 
measured to the computed g^ from the k-iterated solution, and 

is the normalized weighting function 8P(h,h^)/9h where K.^ <_ 1 . The 
iterative process is terminated when the computed radiant approaches the 
measured radiant to within the instrument noise level. 


Inversion of Multiwavelength Extinction Profiles into Constituent Profiles 


The inversion methods discussed in the previous section will generate 
multispectral extinction vertical profiles from the multichannel extinction 
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measurements. The second step in the inversion process will be to separate 
the contributions from individual constituents. For the wavelength region of 
interest, the total extinction coefficient at u fixed altitude level can be 
written as the sum of contributions £rc\t aerosol, ozone, and air molecules as 

0 total “ ^aerosol ^ + a Ai P mol C V + b Ai 0 0 3 ( V (13) 

whore S mol (A Q ) is the Rayleigh extinction coefficient at a reference wave- 
length A , a. . is the wavelength dependence of the Rayleigh extinction, 
b^ is the wavelength dependence of the ozone adsorption profile, and 
^aerosol ^ tbe aeroso1 extinction as a function of wavelength A , 

The Rayleigh extinction versus wavelength profile is given by the fourth 
power law a.^ = (A Q /A^) ,t , The ozone absorption profile for the Chappuis 
band has been measured (Inn and Tanaka, IPS?), The aerosol extinction versus 

wavelength dependence is not known and is dependent on aerosol parameters 

* 

such as shape, size distribution, and complex index refraction. If we 
assumed the aerosols are all homogeneous spheres, then for a given size 
distribution N(r) , the aerosol extinction profile is given by 


Aerosol CX) = / CT(X ’ r} N(r} dr Cl4) 

where o(A,r) is the scattering cross section for the spherical aerosol 
of radius r at wavelength A . The inversion of equation (14) to solve for 
the size distribution N(r) , even assuming the complex index of refraction 
is known, is a very difficult task (Twomey and Howell, 1967). In this paper, 
we will not attempt to do this inversion. Instead, we will assume that a 
two-parameter model will be sufficient for the description of aerosol wave- 
length dependent extinction coefficient. We assume the two-parameter model 
to be of the following form 

Aerosol ' W * A t 15 > 
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where A and a are the two paramours. The functional form of the aerosol 
extinction coefficient as expressed by equation (15) is identical to the result 
obtained by assuming the aerosols to have a Junge size distribution (Junge, 
1963). Substituting equation (15) into (13), we then have a system of nonlinear 
equations with unknown 8 mo j.(A Q ) > ^03^0^ ‘ ^ * anc * a ^ or eac ^ a ^ t: ^ tu ^ 0 

level. The system of nonlinear equations can be solved using the Marquardt 
(1963) algorithm. The Marquardt algorithm is a minimum search procedure for a 
system of nonlinear equations. An initial set of guessed solutions for the 
unknowns is updated through the following iterative equation 

X (k+D . x (k) _ [j T ( x k ) j( x k ) + x kI ]-l j T ( x k ) p(x k ) (16) 


k 

where X is the solution vector for the unknowns at the k-iteration, 

F(x k ) is the residue vector for equation (13) at the k-iteration, and 
J(x k ) is che Jacobian matrix with element J.. = 3f./3x. where f. is the 
ith nonlinear equation. I is the identity matrix, and is a control 

parameter at the k-iteration. The control parameter is adjusted at each 
iteration for fast convergence to the minimum. The iterative process is 
terminated when the difference between two consecutive residue squares is 
less than 0.1 percent. 


INVERSION RESULTS 


Inversion results have been obtained by applying the inversion techniques 
discussed above to simulated spacecraft solar extinction measurements. 

Several models of the aerosol and ozone profiles were used to generate 
atmospheric transmittance profiles at different wavelength regions. Wavelengths 
at 0.38, 0.45, 0.6, and 1.0 micron were chosen to coincide with the four 
spectral channel locations for the Stratospheric Aerosol and Gas Experiment 
(SAGE) . The simulated radiant intensity data are generated from the 
atmospheric transmittance profiles by incorporating the spacecraft geometry. 

The radiometer is assumed to have a field of view of 30 arc seconds situated 
on a satellite at orbital altitude of 600 km. A total of 80 radiant intensity 
data points are generated per channel, covering tangent altitudes from 10 to 
50 km in equal increments of half km height. Each data point is simulated by 
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assuming the radiomotor is pointing to a different location on the solar 
disk. 

The inversion of these radiant profiles was then performed by intro- 
ducing experimental errors of various origin and magnitudes to the simulated 
data. There are three types of errors which directly affect the inversion 
results. The first type of error is the bias error in the determination of 
the exact tangent heights for all the data points. The second type of error 
is the noise associated with each measurement. The third type of error is 
associated with the uncertainty in knowing the pointing angle of the radio- 
meter for each data point. This pointing error will result in both an 
erroneous value for the unattenuated solar radiants due to solar limb 
darkening effect and an erroneous value for the instantaneous tangent height. 

The three types of errors mentioned above in different magnitudes were 
added to the simulated radiant data, and inversions were performed to deter- 
mine the sensitivity of these errors on the inversion results. The atmosphere 
is divided into 45 hoi 5 .v w en^ous layers of half km thickness between 10 and 21 
km, and increasing thickness from 21 to 50 km. The path length matrix is 
computed using a ray trace calculation to include atmospheric refraction 
effect. 

Figures 3 to 6 show inversion results for simulated extinction measure- 
ments using the linear constrained inversion method. The input aerosol 
vertical profile used is deduced from an actual lidar backscattered profile 
with vertical resolution of 0.15 km. The ozone vertical profile used is 
similar to Elterman's model (1968) except that artificial layers structures 
were added. The Rayleigh extinction profile is also deduced from an actual 
rawinsonde temperature and pressure profile. In these simulated 
measurements experimental errors including noise level of 0.1 percent to full 
scale radiant intensity and pointing uncertainty of 3 arc seconds were added 
to the simulated measurements. The inversion results are presented in terms 
of aerosol vertical extinction at 1.0 micron (figure 3), ozone vertical 
extinction at 0.6 micron (figure 4), Rayleigh vertical extinction at 0.38 
micron (figure 5), and aerosol optical model parameters (figure 6). The 
aerosol optical model parameter a is assumed to be -0.8 between 17 and 
20 km and equal to -1.6 for the rest of the altitude range. Inspection of 
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Figure 3. Inversion results for aerosol vertical extinction profile at 1.00 micron 
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Figure 6. Inversion results for aerosol optical model parameter a as a function of 
altitudes. 



figures 3 through 6 indicates that both aerosol and ozone vertical profiles 
can be inverted with good accuracy up to 40 km altitude. The poor inversion 
results for the Rayleigh profile below 14 km are caused by the rapidly decreasing 
signal level at the 0.38 micron channel. In figure 6, the inverted results for 
the aerosol optical parameter a within the Junge layer region show scm, 
fluctuation, and no inversion results can be obtained above 20 km due to the 
rapid decrease of aerosol concentration in this region. 

Figures 7 to 9 show inversion results for the same input models except 
that the noise level for each measurement has been increased to 1.0 percent, 
and the pointing uncertainty has been increased to IS arc seconds. Both the 
inverted results for aerosol and ozone profiles show considerable decrease 
in resolution. This behavior is consistent with the inversion technique in 
which the smoothing parameter a Q in equation (11) is increased to 
accommodate the high noise level in the measurements. The inversion results 
for the aerosol optical model parameter which is not shown here are similar to 
figure 6 except larger fluctuations were observed. 

Figures 10 and 11 show inversion results for different input aerosol and 
ozone vertical profiles. The aerosol profile is deduced from the first 
published lidar observation of the Volcano de Fuego eruption (McCormick et al., 
1975) showing sharp layering structures. The peak layer at 19 km is approxi- 
mately 0.5 km wide. The inversion cannot reproduce this sharp layering 
structure due to the integrating effect over the radiometer's field of view. 

Figures 12 and 13 show the inversion results using the iterative method 
as discussed in equation (12) . The atmospheric model and experimental noise 
levels are identical to those shown in figures 3 to 6. Comparison of the 
inversion results in this case to those shown in figures 3 and 4 indicates 
that they are quite similar. The inverted aerosol profile at altitudes above 
30 km in this case shows some large amplitude oscillations. This is expected 
as the noise level in this altitude region is considerably higher than the 
signal level. 

Inversions have also been performed on simulated data including bias 
errors in the tangent height determination. It is found that bias errors will 
produce a shift in altitude scale of the inverted profile in direct proportion 
to the bias magnitude. 
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Figure II. Sane as figure 4 with different ozone vertical extinction profile 
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CONCLUSIONS 


This report has demonstrated that measurements from a spacecraft solar 
extinction experiment can be inverted to produce aerosol and ozone vertical 
profiles from cloud top up to approximately 50 kin altitude. Analysis of the 
inversion results from simulated measurements including various experimental 
errors indicated that the resolution of the inverted profiles will be 
degraded as the errors are increased. 

Both the linear constrained inversion method and the iterative method 
have been used for the inversion of solar extinction measurements. The 
accuracy of the inverted results from the two different inversion methods 
are similar. 
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